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ABSTRACT 

The effects of strong gravity on the polarization of the Compton reflection from an 
X-ray illuminated accretion disc are studied. The gravitational field of a rotating black 
hole influences Stokes parameters of the radiation along the propagation to a distant 
observer. Assuming the lamp-post model, the degree and the angle of polarization are 
examined as functions of the observer's inclination angle, of the height of the primary 
source and of the inner radius of the disc emitting region. It is shown that polarimetry 
can provide essential information on the properties of black holes sources, and it is 
argued that time variation of the polarization angle is a strong signature of general 
relativity effects. The expected polarization degree and angle should be detectable by 
new generation polarimeters, like the one planned for the Xeus mission. 
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1 INTRODUCTION 

Accretion discs in central regions of active galactic nuclei 
(AGN) are subject to strong external illumination originat- 
ing from some kind of corona and giving rise to specific 
spectral features in the X-ray band, namely a Compton Re- 
flection component and fluorescent emission lines, the most 
prominent of them being that of iron. It has been shown 
that the shape of the intrinsic spectra must be further mod- 
ified by the strong gravitational field of the central mass, 
and so X-ray spectroscopy could allow us to explore the in- 
nermost regions of accretion flows near super-massive black 
holes (for recent reviews, see Fabian et al. 2000; Reynolds 
& Nowak 2003). Similar mechanisms operate also in some 
Galactic black hole candidates. 

However, a rather surprising result from recent XMM- 
Newton observation is that relativistic iron lines are not as 
common as previously believed (Bianchi et al. 2004 and ref- 
erences therein). This does not necessarily mean that the 
iron line is not produced in the innermost regions of ac- 
cretion discs. The situation is likely more complex than in 
simple, steady scenarios, and indeed some evidence for the 
line emission arising from orbiting spots is present in the 
time-resolved spectra of a few AGN (e.g., Dovciak et al. 
2004a, and references therein). Even when clearly observed, 
relativistic lines behave differently than expected. The best 
example is the puzzling lack of correlation between line and 
continuum emission in MCG-6-30-15 (Fabian et al. 2002), 
unexpected because the very broad line profile clearly indi- 



cates that the line originates in the innermost regions of the 
accretion disc, hence very close to the illuminating source. 
Miniutti et al. (2003, 2004) have proposed a solution to this 
problem in terms of the source moving along the black hole 
rotation axis or very close to it. 

In this paper we demonstrate that polarimetric studies 
could be very useful to discriminate between different ge- 
ometries and physical states of accreting sources in strong 
gravity regime. The idea of using polarimetry to gain addi- 
tional information about compact objects is not a new one. 
In this context it was proposed by Rees (1975) that polarized 
X-rays are of high relevance. Pozdnyakov, Sobol & Sunyaev 
(1979) studied spectral profiles of iron X-ray lines that re- 
sult from multiple Compton scattering. Laor et al. (1990) 
studied the polarization of a thin accretion disc around a 
rotating black hole within the framework of general rela- 
tivity, using an improved radiative transfer computations. 
Later on, various influences affecting polarization (due to 
magnetic fields, absorption as well as strong gravity) were 
examined for black hole accretion discs (Agol & Blaes 1996; 
Ogura, Ohuo & Kojima 2000). Temporal variations of polar- 
ization were also discussed, in particular the case of orbiting 
spots near a black hole (Connors, Piran & Stark 1980; Bao, 
Wiita & Hadrava 1996). Here we examine consequences of 
a specific model of an illuminated accretion disc. With the 
promise of new polarimetric detectors (Costa et al. 2001), 
quantitative examination of specific models becomes timely. 

Since the reflecting medium has a disc-like geometry, a 
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substantial amount of linear polarization is expected in the 
resulting spectrum because of Compton scattering. Polar- 
ization properties of the disc emission are modified by the 
photon propagation in gravitational field, providing addi- 
tional information on its structure. In this paper we calculate 
the observed polarization of the reflected radiation assum- 
ing the lamp-post model (Martocchia & Matt 1996; Petrucci 
& Henri 1997; Ghisellini et al. 2004). We assume a rotating 
(Kerr) black hole as the only source of the gravitational field, 
having a common symmetry axis with an accretion disc. 

In the next section, by employing Monte-Carlo com- 
putations (Matt, Perola & Piro 1991; Matt 1993) we will 
find the intrinsic emissivity of an illuminated disc. Then we 
will integrate contributions to the total signal across the 
disc emitting region and, using a general relativistic ray- 
tracing technique (Dovciak et al. 2004b, 2004c), we derive 
the Stokes parameters as measured by a distant observer. We 
will present the polarization properties of scattered light as 
a function of model parameters, namely, the height z = h 
of the primary source on the symmetry axis, the dimension- 
less angular momentum a of the black hole, and the viewing 
angle o of the observer. The conclusion will be that the 
X-ray polarimetry can be very useful to help discriminating 
between models and determining their parameters. 



2 POLARIZATION OF REFLECTED LIGHT IN 
THE STRONG GRAVITY REGIME 

2.1 An irradiated disc near a black hole 

The local emissivity of radiation reflected at a given point of 
the disc is assumed to be proportional to the incident illumi- 
nation by the primary power-law continuum. The primary 
source of the lamp-post model is located at height h above 
the black hole, and so we first need to integrate light rays 
from the source down to the disc. The disc is assumed to 
be stationary and we restrict ourselves to the time-averaged 
analysis, assuming processes that vary in a much slower pace 
than the light-crossing time at the corresponding radius. 

An incident photon strikes the disc at a certain point 
V(r,ip) in equatorial plane 8 = 7r/2 (i.e. z — 0). Four- 
momentum of an incident photon is 
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) A p — a is Carter's constant 
of motion, A(r)=r 2 - 2r + a 2 , A p =A(h), and TL a is the 
sign of p\. Usual notation is employed for Kerr metric func- 
tions in Boyer-Lindquist coordinates (e.g. Kato, Fukue & 
Mineshige 1998). All quantities are made dimensionless by 
scaling with gravitational radius GM/c 2 . The disc is defined 
by r- m <r<r out . We also denoted 8 P to be the local angle 
of emission under which photon emerges in the source rest 
frame (9 P — corresponds to a photon heading downwards 
to the disc, while 6 P = 180° is for upward direction). Affine 
parametrization has been employed in such a way that con- 
served energy of an incident photon is — pi t = — p p t = 1 and 
its conserved angular momentum vanishes, i.e. l p = 0. The 



assumed range of key parameters is 0<a<l, 0<# o <90° (in 
this paper we show results for a rapidly rotating black hole, 
a — ► 1, and we use 5° resolution in o ). 

The following considerations are necessary in order to 
derive the observed spectrum and polarization. First, the 
gravitational and Doppler-induced shift of energy of the pri- 
mary photons impinging on the disc is 
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where E p and E\ denote the photon energy at the point 
of emission from the primary source and at the point of 
incidence on the disc, respectively. (7 M is Keplerian four- 
velocity of the disc medium and U p is four- velocity of the 
primary source. The only non-zero component of the latter 
quantity is U l p = ^ (h 2 + a 2 ) A p 1 . 
Cosine of local incidence angle is 
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, 0) denotes normal direction 



to the disc in the co-moving frame of the disc medium. For 
explicit definition of the local tetrad ef a \ and further details, 
see Dovciak et al. (2004b). 

The reflection component has been computed by a 
Monte-Carlo code (Matt 1993). The number of reflected 
photons is proportional to the incident flux N-f (E p ) arriving 
from the primary source, 



N f{E p ) = N?(E p )£±- 

ClOloc 



(7) 



where iVp (E p ) — No p E p T represents an isotropic and 
steady power-law primary emission that is emitted into solid 
angle df2 p and eventually illuminates the local area element 
dSioc on the disc. The ratio df2 p /d£'i oc is 
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where dS = r~ 2 (pf)~ 1 dS± is the element of coordinate area 
corresponding to proper area dS± perpendicular to the in- 
cident ray. The proper area projected onto the disc is 
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Only direct photon rays are considered in the present com- 
putations, indirect image photons being known to pro- 
duce merely a marginal correction to the total signal under 
usual circumstances, i.e., a moderate inclination angle of the 
source. This is because luminosity of the successive images 
of the same source decreases exponentially with the order of 
the corresponding image (e.g., Luminet 1979). The contri- 
bution of indirect photons may be important only in case of 
a highly inclined (edge-on) disc system, in which case caus- 
tics play a role (Bao et al. 1994; Rauch & Blandford 1994). 
It follows from equations that the incident spectrum 

on the disc conforms to a power-law profile in energy with 
the same photon index F as the original primary emission, 
i.e. Nf (E) = Aoi E~ r with the normalization factor 



A 0i = N 0p g p 
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Figure 1. Polarization degree and angle due to reflected radiation integrated over the whole surface of the disc and propagated to 
the point of observation. The dependence on height h is plotted for three values of observer's view angle O (edge-on view of the disc 
corresponds to O = 90°). Left panel: r; n = 6; right panel: r; n = 1.20. In both panels the energy range was assumed 9-12 keV, photon 
index of the incident radiation T = 2, angular momentum a = 0.999. 
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Figure 2. Polarization degree and angle as functions of fi (cosine of observer inclination) for the same model in previous figure. 



2.2 Gravity induced changes of polarization 

The assumption of a geometrically thin accretion disc near a 
rotating black hole defines geometry of the system and met- 
ric of the space-time. The intrinsic polarization of emerging 
light can be then computed locally, assuming a plane-parallel 
scattering layer. This problem was studied extensively in 
various approximations (Chandrasekhar 1960; Sunyaev & 



Titarchuk 1985). Here we adopt the approach of Matt (1993) 
and we refer the reader to recent papers (Dovfiiak et al. 
2004b, c) for more details on the treatment of polarized light 
in strong gravitational fields. 

Four Stokes parameters, /„, Q v , U v and V„, entirely 
describe polarization properties of the scattered light. Here- 
after we will distinguish the quantities that are determined 
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locally at the point of emission on the disc surface (denoted 
by index 'loc') and those relevant to a distant observer (in- 
dex 'o'). We introduce specific Stokes parameters, 



E 



Uu 

E 



Vv 
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(11) 



and then specific Stokes parameters per energy bin, i.e. Ai , 
Aq a , Aw and Av a . The latter quantities are directly mea- 
surable, specifying the fluxes of photons with a given polar- 
ization. We can write 

Ai (E,AE) = No J dS J dEi oc ii oc (Ei oc )F , (12) 

Aq {E,AE) = No J dS J dE ioc [q loc (E loc ) cos 2* 

-wioc(-B loc )sm2*] J F I (13) 

Au (E,AE) = No J dS J dE ioc [gioc(#loc) sin2* 

+Mioc(Si oc )cos2*]F, (14) 

Av (E,AE) = No dS dE loc v loc {E loc )F. (15) 



Here, F = F(r,tp) = g 2 lfi c r is the transfer function, g 
being the total energy shift between observed and emitted 
photons, I the lensing effect, /x e the cosine of the emission 
angle, and \& the angle by which a vector rotates while it is 
parallelly transported along the light geodesic. One can refer 
to $ as the gravity- induced change of polarization angle, be- 
cause polarization vector is parallelly transported along the 
light ray. The integration domain covers the X-rays emitting 
surface of the disc, ri n <r<r ou t- Notice that the local specific 
Stokes parameters may depend on r, tp and fi c , which we do 
not state explicitly in eqs. 1121 - 1151 for simplicity but we 
take this dependency into account when integrating contri- 
butions to the observed values. 

A complementary way of characterizing polarization is 
in terms of degree of polarization P and by two polarization 
angles Xo and £ D : 



Po 

tan 2\ a 
sin 2£ 



= \/ ql + u% + vi /io , 
= Uo/q , 

= Vo/y/qi + ul + vi. 
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Unpolarized primary radiation was assumed and the approx- 
imation of single Rayleigh scattering was adopted in our ex- 
amples hereby. The angular dependence of the local Stokes 
parameters is therefore given by eqs. (1.147) and (X.172) 
in Chandrasekhar (1960). See Matt (1993) and Dovciak et 
al. (2004b) for discussion of the adopted approximation in 
the present context. The angle £ characterizes the circular 
polarization and it vanishes in our case. 

An explicit formula can be given for the angle $i be- 
tween the projection of three-momentum of an incident pho- 
ton (in the local rest-frame co-moving with the disc) and the 
radial vector, 



$i = -TV a 



-1 



(19) 



where TV a is positive for incident photons travelling outwards 
(Pi (r) > 0) and it is negative in case of inward direction 



(p; r ' < 0). The (^-component of the local tetrad is e 



(f) a 



>/S(-^,0,0,J7'). Equation JUJ appears in evaluation of 
local Stokes parameters (Chandrasekhar 1960). 

Polarization of scattered light is shown in Figure Q 
where we plot the polarization degree and the change of po- 
larization angle as functions of h (we used a 14 points grid 
with h < 100). Notice that in the Newtonian case only polar- 
ization angles of 0° or 90° would be expected for symmetry 
reasons. The two panels in the figure correspond to differ- 
ent locations of the inner disc edge: ri n = 6 and r ln = 1.20, 
respectively. The curves are strongly sensitive to n n and h, 
while the dependence on r out is weak for a large disc (here 
r ut = 400). Sensitivity to n n is particularly appealing if one 
remembers practical difficulties in estimating r- m by fitting 
spectra. The effect is clearly visible even for h ~ 20. Graphs 
corresponding to r- ln = 6 and a — 0.999, resemble, in essence 
quite closely, the non-rotating case (a = 0) because dragging 
effects are most prominent near horizon. 

Figure shows the polarization degree and angle as 
functions of the observer's inclination. Again, by comparing 
the two cases of different r ln one can appreciate how sen- 
sitive polarization is to details of the flow near the inner 
disc boundary. Furthermore, it is worth noticing that the 
polarization degree has a local maximum for moderate in- 
clination angles (fi a ~ 0.9). This peak occurs mainly due to 
interplay between the intrinsic polarization on the disc and 
special-relativistic aberration of the outgoing photons, while 
general relativistic effects are shaping details of the profile. 

The above-described polarization features are represen- 
tative of the scattering mechanism and the gravitational 
field structure acting on reflected photons. However, in or- 
der to compute directly observable characteristics one has to 
combine the primary continuum with the reflected compo- 
nent. Polarization degree of the resulting signal depends on 
mutual proportion of the two components and also on the 
energy range of observation. In Figure [3] we assumed that 
the irradiating source emits isotropically and its light is af- 
fected only by gravitational redshift and lensing, according 
to the source location at z = h on axis. This results in di- 
lution of primary light by factor ~<?hC l i &>) lh{h, O ), where 
gh = yl — 2h/(a 2 + h 2 ) is the redshift of primary photons 
reaching directly the observer. The whole term is a com- 
bination of special-relativistic transformation and general- 
relativity light bending (factor lh is order of unity and can 
be safely ignored here; cf. Dovciak et al. 2004b). Anisotropy 
of primary radiation may further attenuate or amplify the 
polarization degree of the final signal, while the polariza- 
tion angle is rather independent of this influence as long as 
primary light is itself unpolarized. 



3 CONCLUSIONS 

We examined the polarimetric properties of X-ray illumi- 
nated accretion discs in the lamp-post model. We found that 
observed values of Stokes parameters are expected to be 
rather sensitive of the model parameters. The adopted ap- 
proach provides additional information with respect to tra- 
ditional spectroscopy, and so it has a great potential for dis- 
criminating between different models. It offers an improved 
way of measuring rotation of the black hole, because the ra- 
diation properties of the inner disc region most likely reflect 
the value of the black hole angular momentum. Our calcula- 
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Figure 3. Net polarization degree of the total (primary plus reflected) signal as a function of height h (left panel) and a function of /i 
(right panel). In each panel, curves are parametrized by the corresponding energy range. 



tions are complementary to those of Laor et al. (1990), where 
the general relativity induced change of polarization angle 
and degree were computed for the disc thermal emission. 

While our calculations have been performed assuming 
a stationary situation, in reality it is likely that the height 
of the illuminating source changes with time, and indeed 
such variations have been invoked by Miniutti et al. (2003) 
to explain the primary and reflected variability patterns of 
MCG-6-30-15. Complete time-resolved analysis (including 
all consequences of the light travel time in curved space- 
time) is beyond the scope of this paper and we defer it to 
future work, assuming here that the primary source varies 
on a time-scale longer than the light- crossing time in the 
system. This is a well-substantiated assumption also from 
a practical point of view, since feasible techniques will any- 
way require sufficient integration time (i.e. of the order of 
several ksec). Here, suffice it to note that a variation of h 
implies a variation of the observed polarization angle of the 
reflected radiation. Although the change of polarization an- 
gle and the degree of polarization are sensitive to general 
relativity effects, to certain degree they depend also on other 
assumptions about the source geometry and local physics of 
the source. However, it should be obvious that time-resolved 
polarimetry can set constraints on the models that are sub- 
stantially more stringent than what can be achieved by pure 
spectroscopy. 

New generation photoelectric polarimeters (Costa et al. 
2001) in the focal plane of large area optics (such as those 
foreseen for Xeus) can probe polarization degree of the or- 
der of one percent in bright AGN, making polarimetry, along 
with timing and spectroscopy, a tool for exploring the prop- 
erties of the accretion flows in the vicinity of black holes. 
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